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Abstract 
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temperature is presented within the framework of hquid drop model, including 

the essential finite size effects. Strong parameter dependences of the properties 

are found and discussed. 
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I. INTRODUCTION 



It was argued by Witten that strange quark matter might be the ground state of 
normal nuclear matter at zero temperature and zero pressure, which was later supported by 
the studies of Farhi and Jaffe |0] based on MIT bag model. The existence of stable strange 
quark matter would have some remarkable consequences in cosmology and astrophysics. 
For instance, it has been suggested that some of the dark matter in the Universe could 
possibly exist in the form of strange quark matter which was produced when the Universe 
underwent the quark to hadron transition. The possible transition of a neutron star into 
a strange star was discussed as well In addition, there have been proposals that the 

strange stars might be the sources of gamma rays with very high intensities PJ^. Recently 
Greiner and his collaborators have proposed that small lumps of strange quark matter 
( "strangelets" ) might be produced in ultrarelativistic heavy-ion collisions, and they could 
serve as an unambiguous signature for the formation of quark-gluon plasma. In fact, several 
heavy- ion collision experiments at Brookhaven and CERN are searching for strangelets 0]. 

Most strange quark matter investigations P,PJT0|1 have been performed based on shell 
model or liquid drop model. While calculations in shell model are rather tedious, recently 
Madsen has pointed out that the general structure of shell model results can be ob- 
tained more readily from liquid drop model. On the other hand, strange quark matter at 
finite temperature has attracted some special interests [ll2| , p!3| , since in some cases, such 
as ultrarelativistic heavy-ion collisions, strange quark matter is expected to exist in a hot 
environment. However, these studies of strange quark matter at finite temperature have 
not completed the treatment of finite size effects, which are expected to play a substantial 
role in strangelets possibly appearing in the scope of ultrarelativistic heavy-ion collision ex- 



periments [|14|. The present work is an attempt to apply liquid drop model to strangelets 
at finite temperature, and could be a step towards the studies of phase structure and evo- 
lution of strange quark matter. We will include the important finite size effects (volume, 
surface and curvature contributions) to study the overall properties of strangelets at finite 
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temperature, which might give some clues to the search of strangelets. 

In the hterature P, p!T] , p!5| there have been discussions on parameter dependences of the 
stabihty of strange quark matter at zero temperature. However, it is interesting as well 
to study how the parameter "windows" for stable strange quark matter would become, if 
strange quark matter at finite temperature, e.g. strange stars with a temperature up to tens 
of MeV, is assumed to exist stably in nature. In this work we shall investigate the properties 
of strange quark matter at finite temperature for a wide range of parameters. 



II. FORMALISM 



We consider the strangelet as a gas of up, down, strange quarks, their antiquarks, and 



gluons confined in an MIT bag. In liquid drop model |TT| , the grand potential of the system 
is given by 

n = J2^i + BV. (1) 

i 

Here B is the bag constant, V is the bag volume, and the grand potential of particle species 
i can be written as 



foo (jM. 

TT dfc— iln[l±exp(-(e,(A;)-/i,)/T)], (2) 

JO CLK 



where the upper sign is for fermions, the lower for bosons, fi and T are chemical potential and 
temperature, k and are particle momentum and energy. In a multiple reflection expansion 
p!6| the density of states is given by 



For spherical strangelets V = |vri?^ is the volume of the bag, 5* = AnR"^ is the surface area, 
and C = SnR is the extrinsic curvature of the bag surface. The factor Qi is the statistical 
weight (6 for quarks and antiquarks, and 16 for gluons). 

The surface term for quarks was given by Berger and Jaffe as 
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In particular, the surface term for massless quarks and gluons is zero. It has been shown by 
Madsen that the following ansatz works for the curvature term of massive quarks: 



For gluons the curvature term is [|T8 



After the construction of the grand potential as above, we can readily obtain the thermo- 
dynamical quantities of the system as follows. The net number of quarks, i.e. the number 
of quarks minus the number of antiquarks, can be derived from 

and 



^ = ^E^» (8) 

gives the total baryon number of the strangelet. The free energy F and internal energy E 
of the strangelet are given by 

F = J2{Q, + N,fi,) + BV, (9) 

i 

and 

E = F + TS (10) 

with the entropy S = - (fp)^^ • 

III. NUMERICAL RESULTS 



We can study the ground state properties of strangelets at finite temperature T, i.e. the 
lowest mass state for a given baryon number A and temperature T, by minimizing the free 



energy F with respect to the net number of quarks Nq, Ns and the volume of the system 
V , at fixed baryon number A. However, for the sake of illustration, in Fig. 1 we vary the 
strangeness fraction fs = Ng/A (net number of strange quarks per baryon) instead of A^^, 
to show the minimum of the free energy per baryon. While the free energy per baryon 
decreases as temperature rises up from zero to 30 MeV, the value of fs which minimizes the 
free energy per baryon changes little. Actually, as it will be shown later, the strangeness 
fraction of strangelet in its ground state has weak dependence on temperature. 

An easier way to investigate the ground state properties of strangelet is to minimize the 
free energy analytically with respect to Ng, Ng and V, under the constraint Eq. (§) for fixed 
A, which leads to the mechanical equilibrium condition 

and the optimal composition condition 

Hg = Us (12) 

(For details see Appendix). Eq. (^) has the explanation that the pressure exerted by 
the quarks and gluons in the bag is balanced by the bag pressure B. Eq. (|I^) might be 
unexpected, however. It could be interpreted as follows. Since the total number of quarks 
(= Ng + Ns = 3 A) in the strangelet is fixed, the free energy of the system is minimized 
when adding a massless quark (up or down quark) to the strangelet is energetically as 
favorable as adding a massive strange quark, i.e. /i^ = /i^. Eq. (|12D is usually considered as 
a consequence of bulk strange matter at energy minimum. We argue that Eq. ([T2|) holds for 
strangelets as long as the liquid drop model is applicable (applicability of liquid drop model to 
strangelets has been approved by Madsen |11|). Nevertheless, if some other contributions, 
e.g. Coulomb effects, are taken into account, the simple relationship in Eq. (|1^) will be 
modified. Eqs. (H) — (|T2D altogether specify the ground state properties of strangelets. 

A strangelet at zero temperature is stable relative to ^^Fe nuclei if its energy par baryon, 
E/A, is less than 930 MeV. At finite temperature, however, a strangelet can radiate ele- 
mentary particles like nucleons and pions from its surface. While it does not lead to the 



breakup of a hot strangelet but rather favors the creation of a cold strangelet, pion emis- 
sion is suppressed at low temperature by a factor ~ exp{—mT^/T) (~ 0.01 at T=30 MeV). 
Furthermore, for a strangelet with E/A <930 MeV, nucleon emission proceeds via energy 
fluctuations and will be greatly suppressed ||7|,p!2|. If E/A >930 MeV, on the other hand, 



nucleon emission process has no threshold and goes much faster, which might drive the hot 
strangelet to a complete hadronization. Therefore, in this context we shall call a strangelet 
stable if its E/A is lower than 930 MeV throughout this work. 

In Fig. 2(a) we have presented the (internal) energy per baryon as a function of tem- 
perature, for two different baryon numbers, A = 20 and 100, with bag constant 5^/^=145 
MeV and strange quark mass ms=150 MeV. Note that the strangelet with a baryon number 
A=100 is stable only for temperature T <20 MeV. (The dot-dashed line in Fig. 2(a) marks 
an energy per baryon of 930 MeV). We see also that the upper limit of temperature for 
stable configuration of strangelets depends on baryon number. This is a consequence of the 
inclusion of finite size effects. More discussions will be given in Fig. 4. 

Fig. 2(b)-(d) show the temperature dependence of the strangeness fraction fg, charge-to- 
baryon-number ratio Z/A, and averaged radius per baryon R/A^^^ of strangelets with A=20, 
and 100, for 5^/^=145 MeV and 771^=150 MeV. It can be seen that fs and R/A^^^ increase, 
and Z/A decrease slowly, when the temperature rises up. It is expected that fs will saturate 
to unity, and Z/A will decrease to zero, when the temperature is high enough, since at that 
stage the mass difference between up, down and strange quarks is unimportant, the system 
will tend to become flavor symmetric (A^^ = = Ng) and globally charged neutral. Our 
observation of Z/A as a decreasing function of temperature differs from that of Chakrabarty 

In Fig. 3(a)-(d) we have plotted the energy per baryon E/A, strangeness fraction /<,, 
charge-to-baryon-number ratio Z/A, and averaged radius per baryon R/A^^^ as a function 
of baryon number A, for two temperatures T=0 and 30 MeV, with 5^/^=145 MeV and 
m<j=150 MeV as well. One can see from Fig. 1(a) that there is a lower limit A^ri to the 
baryon number of a stable strangelet (whose E/A <930 MeV). This critical baryon number 
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is a function of temperature, parameters B (bag constant) and (strange quark mass), 
as will be shown later. Except for sufficiently low baryon number, we see that fs, Z/A and 
R/A^/^ depends weakly on A. Their typical values are fs ~0.4-0.6, Z/A ~0. 15-0.3, R/A^^^ ~ 
0.95-1 fm. 

Fig. 4 shows the critical baryon number as a function of temperature, for two different 
values of mg, =150 MeV (full line) and ms=0 (dashed hne) with 5^4=145 MeV. For 
temperatures low enough, the critical baryon number increases little. However, beyond 
some temperature (about 30 MeV for ms=150 MeV, and 42 MeV for nis—O) even bulk 
strange quark matter (^4 —>■ oo) could not be stable. 

Before any concrete statements can be made, it is reasonable to discuss the parameter 
dependence. Next we study how the properties of strangelets depend on the two main 
parameters in liquid drop model, bag constant B and strange quark mass m^. 

The variation with S^/^ are plotted in Fig. 5(a)-(d) for strangelets with a baryon number 
^4=50, at temperature T—0 (full lines) and 30 MeV (dashed lines), with ms=150 MeV. One 
can see from Fig. 5(a) that the energy per baryon E/A is proportional to B^^^ for a given 
baryon number and temperature. While the strangeness fraction fg increases with 5^/^, the 
charge-to-baryon-number ratio Z/A and averaged radius per baryon R/A^^^ are decreasing 
functions of i?^/^. 

In Fig. 6 the critical baryon number Ach is plotted as a function of S^/^, with the strange 
quark mass 777,5=150 MeV. It can be seen that for sufficiently large bag constant, there could 
be no stable strange quark matter. This upper limit of B^^^ depends on temperature, 
<153 MeV at r=0, < 151 MeV at r=15 MeV, and <145 MeV at 7=30 
MeV. For strange quark mass less than 150 MeV, the upper limit of B^/^ moves up. 

Fig. 7(a)- (d) show the variation with of properties of strangelets at temperature T=0 
and 30 MeV, with baryon number A—bO and 5^/^=145 MeV. For strangelets at T—0, both 
the energy per baryon E/A and charge-to-baryon number ratio Z/A increase with increasing 
rUs, and saturate when reaches the chemical potential ~ fis- The strangeness fraction 
fs is a decreasing function of m^, and goes to zero when 777^ ^ /is- For finite temperature, 
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T=30 MeV, the saturation happens at a larger value of m^. However, the averaged radius 
per baryon R/A^/'^ depends little on m^. 

We have plotted in Fig. 8 the critical baryon number Ach as a function of m^, with 
5^/^=145 MeV. Note that at zero temperature Acri increases with increasing m^, and sud- 
denly becomes a constant at ~ /i^, indicating the fact that for ^ /i^, no strange 
quarks exist (see Fig. 7(b) where /s — 0), and an increased value of does not lead to an 
increase of strangelet energy. 



IV. DISCUSSIONS AND CONCLUSIONS 

In the present work we have neglected QCD radiation corrections, as argued in Refs. 
PJT7[| that an increased value of the strong coupling constant could be absorbed into a 



decrease in the bag constant B. Coulomb effects could be included in liquid drop model 



self-consistently as shown by Madsen |]T5[. We show also in Fig. 1 by the dashed lines the 
free energy per baryon versus strangeness fraction after the inclusion of Coulomb energies. 
One sees that Coulomb effects are rather negligible for the strangelet energy, but could 
cause a change of the quark composition of the strangelet (Note in Fig. 1 that the value 
of fs minimizing F/A is shifted a little by the inclusion of Coulomb effects). The resultant 
corrections to Z/ A are within the uncertainties caused by parameters B and m^, however. 

We have found in this work that the stability of strangelets at zero and finite temperature 
depends strongly on parameters 5, the bag constant, and mg, the strange quark mass. An 
increased value of B or nis (for < /i^) tends to destabilize the strangelet. For ms=150 
MeV, the strangelet can be stable only for 5^/^ <153 MeV at T=0, 5^/^ <151 MeV at 
T=15 MeV, and B^'^ <145 MeV at r=30 MeV. For 5^/^=145 MeV, the strangelet is stable 
only for T <30 MeV with ms=150 MeV, and T ^42 MeV with ms=0. This can impose 
important restrictions on the possible existence of hot strange quark matter, e.g. strange 
stars in astrophysics. From the fact that ordinary nuclei are composed of nucleons instead 
of up-down quark matter we can obtain a lower limit [B^/^ ^145 MeV) of the bag constant 
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15| . For the most extreme choice of parameters, i.e. 5^/^=145 MeV and ms=0 (see dashed 
hne in Fig. 4), we find that strange quark matter at a temperature T ^42 MeV cannot exist 
stable. 

The critical baryon number Ach, below which the strangelet cannot be stable, is an 
increasing function of T, B and m^. For ms=150 MeV and 5^/^=145 MeV, the critical 
baryon number is about A^ri —25 at T=0, and A^ri —60 at T=20 MeV. 

The strangeness fraction of the strangelet, fg, is found to increase slowly with the increase 
of temperature T or baryon number A (except for A ^100). However, fg as an increasing 
function of B^^^ and a decreasing function of m^, is sensitive to the values of B^^^ and m^. 
For m,=150 MeV and B^/'^=U5 MeV, the typical values of are about 0.4-0.7. In terms 
of the simple relationship Z/A = (1 — fs)/'^, the charge-to-baryon-number ratio Z/A has 
typical values about 0.15-0.3, for m,=150 MeV and B'^/'^=U5 MeV. 

Finally the averaged radius per baryon, R/A^^'^, is a slowly increasing function of T or 
rris, and a slowly decreasing function of A, whereas it decreases fast with the increase of 
5^/^ For m,=150 MeV, 5^/^=145 MeV, R/A'^/^ ranges from 0.9 to 1.1 fm. 

To conclude, we have shown in the present work that liquid drop model including finite 
size effects provides a successful description of properties of strangelets at zero and finite 
temperature, and allows studies for a wide range of parameters. It is possible to carry out 
further explorations of phase structure, phase evolution, and possible survival of strangelets 
in ultrarelativistic heavy ion collisions within the framework of liquid drop model. 
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APPENDIX: 

Eqs. (|ll]) and (^) are given by minimizing the free energy 
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F=Y.{n, + iV,/i,) + BV, (Al) 

i=q,s 

with respect to V, Ng and A^^, under the constraint 

J2Ni = 3A (A2) 

i 

for fixed baryon number A. For the sake of simphcity we omit the surface and curvature 
terms in F in the following derivations. 
Constructing an auxiliary function 

F'=J2 {^riV,i^,iV,Ni)) + N,i^,{V,N,)} + BV + x( ^ N, - 3a] , (A3) 

i=q,s \i=q,s j 

we can derive the constrained minimum from 



and 

IdF' 



= 0, (A4) 

Ni 



0. (A5) 



or 



V 

Since {dVLi/diii)v = -Ni, Eq. leads to 



EfS) +B = 0. (A7) 



dV 

■fc \ 

which is Eq. (|ll]). 

From Eq. (|A5|) we obtain 



Mi 



or 



/i^ + A = 0, (A9) 

which implies 

yUg = /is- (AlO) 

This is Eq. (|l|). 

10 



REFERENCES 

[1] E. Witten, Phys. Rev. D 30, 272 (1984). 

[2] E. Farhi and R.L. Jaffe, Phys. Rev. D 30, 2379 (1984). 

[3] J. Madsen, Phys. Rev. Lett. 61, 2909 (1988). 

[4] C. Alcock and A. Olinto, Ann. Rev. Nucl. Part. Sci. 38, 161 (1988). 

[5] C. Alcock, E. Farhi and A. Ohnto, Phys. Rev. Lett. 57, 2088 (1986). 

[6] P. Haensel, B. Paczynski and P. Amsterdamski, Astrophys. J. 375, 209 (1991). 

[7] C. Greiner, P. Koch and H. Stocker, Phys. Rev. Lett. 58, 1825 (1987); Phys. Rev. D 
44, 3517 (1991). 

[8] J. Barrette et al. , Phys. Lett. B 252, 550 (1990); M. Aoki et al. , Phys. Rev. Lett. 69, 
2345 (1992); K. Borer et al. , Phys. Rev. Lett. 72, 1415 (1994). 

[9] C. Greiner, D.H. Rischke, H. Stocker and P. Koch, Phys. Rev. D 38, 2797 (1988). 

[10] E.P. Gilson and R.L. Jaffe, Phys. Rev. Lett. 71, 332 (1993). 

[11] J. Madsen, Phys. Rev. D 50, 3328 (1994). 

[12] H. Reinhardt and B.V. Dang, Phys. Lett. B 202, 133 (1988); T. Chmaj and W. Slomin- 
ski, Phys. Rev. D 40, 165 (1989). 

[13] S. Chakrabarty, Phys. Rev. D 48, 1409 (1993). 

[14] J. Madsen, Phys. Rev. Lett. 70, 391 (1993). 

[15] J. Madsen, Phys. Rev. D 47, 5156 (1993). 

[16] R. Balian and C. Bloch, Ann. Phys. 60, 401 (1970); T.H. Hansson and R.L. Jaffe, Ann. 
Phys. 151, 204 (1983). 

[17] M.S. Berger and R.L. Jaffe, Phys. Rev. C 35, 213 (1987); 44, 566(E) (1991). 

11 



[18] R. Balian and C. Bloch, Ann. Phys. 64, 271 (1970). 



12 



FIGURE CAPTIONS 



Fig. 1. Free energy per baryon F/A as a function of the strangeness fraction fg with 
(dashed hnes) and without (full lines) Coulomb corrections, at temperature T=0 (upper two 
curves) and r=30 MeV (lower two curves), for strangelets with a baryon number yl=100, 
with m,=150 MeV and B^'^^Uh MeV. 

Fig. 2.(a)Internal energy per baryon E/A, (b) strangeness fraction /g, (c)charge-to- 
baryon-number ratio Z/A^ and (d)averaged radius per baryon R/A^/^ as a function of 
temperature T, for A—2{) (full lines) and A=100 (dashed lines), with ms=150 MeV and 
SV4=i45 MeV. 

Fig. 3.(a)Internal energy per baryon E/A, (b) strangeness fraction /s, (c)charge-to- 
baryon-number ratio Z/A, and (d)averaged radius per baryon R/A^/^ as a function of baryon 
number A, at temperature T=0 (full lines) and 7=30 MeV (dashed hnes), with ms=150 
MeV and 5^4=145 MeV. 

Fig. 4. Critical baryon number Acri as a function of temperature T, for two different 
values of m^, rus =150 MeV (full line) and ms=0 (dashed hue), with 5^4^145 MeV. 

Fig. 5.(a)Internal energy per baryon E/A, (b) strangeness fraction /s, (c)charge-to- 
baryon-number ratio Z/A, and (d)averaged radius per baryon R/A^/^ as a function of bag 
constant B^/^, for strangelets with baryon number ^4=50, at temperature T=0 (full lines) 
and T=30 MeV (dashed lines), with m,=150 MeV. 

Fig. 6. Critical baryon number Ach as a function of bag constant S^/^, for strangelets at 
temperature T=0 (full hue), r=15 MeV (dot-dashed line), and r=30 MeV (dashed line), 
with m,=150 MeV. 

Fig. 7.(a)Internal energy per baryon E/A, (b) strangeness fraction fs, (c)charge-to- 
baryon-number ratio Z/A, and (d)averaged radius per baryon R/A^^^ as a function of strange 
quark mass m^, for strangelets with a baryon number A=50, at temperature T=0 (full lines) 
and r=30 MeV (dashed hnes), with 5^4=145 MeV. 

Fig. 8. Critical baryon number Acri as a function of strange quark mass m^, for strangelets 
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at temperature 7=0 (full line) T=15 MeV (dot-dashed line) , and 7=30 MeV (dashed line), 
with 5^4=145 MeV. 
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